m 

o 

(N 

Oh: 

< 



Angular dependence and symmetry of Rashba spin torque in ferromagnetic 

heterostructures 

Christian Ortiz Pauyac/ Xuhui Wang,^ Mairbek Chshiev,^ and Aurelien Manchon^'Lj 
^ King Abdullah University of Science and Technology (KAUST), 



Physical Science and Engineering Division, Thuwal 23955-6900, Saudi Arabia 
^SPINTEC, UMR CEA/CNRS/UJF-Grenoble 1/ Grenoble- IN P, INAC, Grenoble, F-38054, 

(Dated: April 18, 2013) 



France 



In a ferromagnetic heterostructure, the interplay between a Rashba spin-orbit coupling and an 
exchange field gives rise to a current-driven spin torque. In a realistic device setup, we investigate the 
Rashba spin torque in the diffusive regime and report two major findings: (i) a nonvanishing torque 
exists at the edges of the device even when the magnetization and effective Rashba field are aligned; 
(ii) anisotropic spin relaxation rates driven by the Rashba spin-orbit coupling assign the spin torque 
a general expression T = T^{e)m x (y x m) -f Tl{e)y x m -f r|(e)m x (z x m) + Tl{9)z x m, 
where the coefficients T^'^ depend on the magnetization direction. Our results agree with recent 
experiments. 

PACS numbers: 75.60.Jk 
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The concept of current-driven spin-orbit torque in 
ultrathin ferromagnetic heterostructuresi*^ and diluted 
magnetic semiconductors^ is attracting much attention 
for providing an efficient magnetization switch mecha- 
nism using just one ferromagnet. In contrast to the 
conventional spin-transfer torque that demands a spin- 
polarized current generated by a reference ferromagnet 
(polarizer) fii§. the spin-orbit torque accomplishes mag- 
netization switching by transferring angular momentum 
between the spin and orbital degrees of freedom through 
a spin-orbit coupling. In ferromagnetic heterostructures 
typically made of magnetic trilayers comprising an ul- 
trathin ferromagnetic film sandwiched between a noble 
metal and an insulator, experiments and theories have 
uncovered a spin-orbit torque of the form^i^ii 



not captured by the earlier models based on either spin- 
Hall effect or Rashba torques f^^ 



Tiim X (y X m) + T^y x m, 



(1) 



where m is the magnetization direction and y is the di- 
rectional unit vector (see Fig. 1). Two components in 
Eq.ilT]) are usually referred to as in-plane (Tii) and per- 
pendicular (Ti) torques. 

The current understanding of the origin of spin-orbit 
torque in ferromagnetic heterostructures combines spin- 
Hall effeclj^ and band structure-induced effects. In the 
former, a spin current is generated by the spin-Hall effect 
in the noble metal layer and injected into the ferromagnet 
to produce a torquci^iiS In the latter, symmetry break- 
ing across the interface between the noble metal and the 
ferromagnet induces a spin splitting in the band struc- 
ture, leading to a nonvanishing current-induced spin- 
orbit field. One version of this band structure-induced 
spin-orbit torque is the so-called Rashba torquei where 
an electric field E embedded across the interface pro- 
duces a nonequilibrium Rashba field Bn (x zxE. Recent 
works have emphasized the complexity of these torques as 
functions of materials parameters.— Of most interest to 
the present study, it has been shown that the spin-orbit 
torque possesses a complex angular dependence that is 
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FIG. 1. (Color online) (a) Schematic view the device cross 
section. A charge current is flowing along the x direction, 
generating an effective Rashba field Br (dotted yellow line) 
that is pointing to y. L is the size of the lateral dimensions, 
(b) Spin torque components T,^ and Te in a spherical coor- 
dinate, with 6 and (p being the azimuthal and the in-plane 
angle, respectively. The coordinates dashed lines show the 
magnetization directions considered in this Letter. 

In this Letter, we employ a set of coupled spin- 
charge diffusion equations^ to characterize the Rashba 
spin torque in a square-shaped device with a width (and 
length) L, see Fig[T] We highlight two major findings. 
First, even when the magnetization m is aligned along 
the effective spin-orbit field (along the y direction), the 
spin-orbit torque does not vanish at the edges of the de- 
vice. Second, the torque amplitudes vary with the mag- 
netization direction (a behavior to be referred as angular 
dependence) and are size-dependent. In the weak (spin- 
orbit) coupling regime, the angular dependence vanishes 
as the device size increases; whereas in a strong coupling 
regime, the angular dependence is insensitive to the de- 
vice size. 

A schematic view of the device is shown in Fig. [Ija). 
The inversion asymmetry across the interfaces gener- 
ates a Rashba type spin-orbit coupling from the poten- 
tial gradient along the z direction. In the quasi-two- 



dimensional system considered here, the diffusive dy- 
namics of nonequihbrium spin density S and charge den- 
sity n are described by^i^^ 

(JTi 

— =DV'^n + BVxy ■ S + TVxy ■ mn + RV^^y ■ m{S ■ m) 



dS .2 s ^ S + SzZ mx{Sxm) 

— — — — jy V ij -— 

Ot Tsf TDP Txc 



(2) 
BV 



xyi 



AxcS X m + 2CVxy X S + 2R{m ■ Vxyn)m 
T[mx {Vxy xS) + Vxy X (m X S)] . (3) 



The parameters arc: fi 
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al^xcVphpT'^/^, R = aA^^rVS, and B = 2a^k%T 
denotes the Rashba spin-orbit coupling and V^^j^ = £ x V. 
The momentum relaxation time is r and spin relaxation 
time due to magnetic impurities is given phenomenologi- 
cally by Tsf., whereas tdp = l/2a^fc|,r is the D'yakonov- 

kpivp) is the Fermi momen- 
= VpT/2 is the diffusion con- 
stant. Briefly, A^c- and C-terms describe spin precession 
around the exchange field and Bji, respectively. B-terni 
couples spin and charge degrees of freedom, leading to 
the electrical spin generation and spin-Hall effect. F-terni 
provides a higher order correction to the preccssional mo- 
tion described by first two terms. i?-tcrm contributes 



Percl relaxation time 
tum (velocity) while D 



to a magnetization renormalization. T^ 



1/AL 



the transverse spin dephasing time in the limit of a 
weak ferromagnet»^ As a boundary condition suggested 
by experimentSf^ the nonequihbrium spin density is re- 
quired to vanish at the transverse edges along the y di- 
rection. 

In a recent work, it was shown that reducing the width 
of the magnetic film dramatically modifies the relative 
magnitude of the in-plane and perpendicular torques in 
the weak Rashba limit.— Here, we further argue that re- 
ducing the size of the device results in changes in the 
symmetry of the torque. To support our argument, we 
plot in Fig. [2] the spatial distribution of the spin torque 
density (along the y axis in the yz— plane) for various 
magnetization directions in both weak (A^c ^ akp) and 
strong {Axe ^ akp) spin-orbit coupling regimes. In Fig. 
[iKb), the spin torque density is expressed in spherical 
coordinates, T = T^e^ + Tgeg, which is more general 
than Eq. ([!]). On the right column of Fig. [51 Tg is ro- 
bust in the bulk, resulting from a robust nonequihbrium 
spin density {Sy) driven by the spin-galvanic effect dis- 
cussed by Edelstein jii This effect disappears towards the 
boundaries, as imposed by the boundary conditions'^ 

An important feature in Fig. [2] appears to be the non- 
vanishing spin torque at the edges even when the mag- 
netization m is parallel to B^ (deep blue curves). In 
general, as the angle between the exchange field m and 
Bn closes, the spin torque amplitude decreases. For a 
strong spin-orbit coupling, the spin-Hall effect drives op- 
positely polarized spin densities Sz accumulating at op- 
posite edgesi^ In our finite-size device, within the dis- 
tance of spin-flip relaxation length from the edges, the 
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FIG. 2. (Color online) Spatial distribution of the Rashba 
torque along the y axis, for a device L = 25 nm. The 
2/2— plane is located at the center of the device, i.e., x = 
12.5 nm. Panels (a) and (c) refer to the in-plane component 
(here Ty = T^). Panels (b) and (d) are for the out-of-plane 
component {T± = Tg). Panels (a) and (b): weak Rashba 
regime (a — 0.001 eV nm, A^c = 0.1 eV). Panels (c) and 
(d): strong coupling (a = 0.05 eV nm, Aj;c = 0.01 eV). The 
inset in panel (c) displays: T^(^rn=y) = jyZL/2'^'t>(^=v)'^vl^ 
for different widths. In panels (a),(b) and (d), T^(^rn=y) is 
multiplied by a factor 10. The Fermi energy is Ep = 0.7 eV, 
kp = 4.3 nm~^ and wf = 5 x 10^'* nm s"^. t = 10"^^ s and 
Tsf = 10-'^ s. 



spin density Sz distributed along the y direction gener- 
ates a nonvanishing local spin torque (at the edges) when 
m is parallel to Bfi. As the spin-orbit coupling weak- 
ens, torques at m II Bpi driven by the spin-Hall effect 
become negligible, see FiglJl^a). Meanwhile, as the sam- 
ple size increases L ^ A^/, the expression T^i^„^^y) = 

{'^/ L) j^~j^,^T^(^m=y)dy exponentially decreases, see in- 
set in Fig. HJc). Another important feature is the inho- 
mogeneous profile of T^ , which is driven by the compet- 
ing effects of spin precession around the total field and 
spin Hall-effect: the spin-Hall effect (precession around 
the total field) is dominant in the strong (weak) Rashba 
regime. When m ^ z, spin density Sz does not con- 
tribute to the spin-Hall induced torque and the in-plane 
torque becomes homogeneous, as confirmed by the solid 
black line in Figl2l(c). In a finite size device, these ef- 
fects contribute to the angular dependence discussed in 
the following. 

First, in an infinite system with a weak Rashba spin- 
orbit coupling, once the anisotropy in spin relaxation 
rates due to the D'yakonov-Perel mechanism is quenched, 
Eq. ^ gives rise to a torque described by Eq. ^^ How- 
ever, as long as the spin relaxation rate is not isotropic, 
the torque assumes a more complex angular dependence. 
In an infinite system, setting '^ xy — z x eEd^x, Eq. (|3]) 
reduces to 

AxcS X m+ T^fn x (S x m) 

-^ xc 

+ — SxX + — SyV + —SzZ = X, (4) 



where the last three terms on the left-hand side subscribe 
to both the D'yakonov-Perel mechanism and spin relax- 
ation induced by random magnetic impurities: T^y = 



'sf 



On the right-hand 



'pp^'^s/ and T^ 1 = 2t^I, 
side, the source term reads 

nE 

X = {By + 2CPy xm + TPmx {y xm)), (5) 

ep 

where P is the spin polarization at Fermi level (ef)- 
Equation (j4]) can be analytically solved in spherical co- 
ordinates using T = T^Brf, -\-Tge.g. In the strong coupling 
limit {B :s> C, r), the spin torque becomes 

T = eg[Py X m + (1 + (^ - /3)^)m xyxm 

ep 

— x(l — oie cos^ 9)[m ■ x)m x z x m 

— aemzinyil + (^ — /?)^)m x z x m 
+ x(^ - /?)(! — ae cos^ 6){m ■ x)z x m 

— l3agmzmyZ X m], (6) 

which comprises one of the major results of this Letter. 
For a succinct discussion, the parameters in Eq.(l6]) are 
defined in Refd. 

Equation ([6]) contains both odd and even components 
with respect to the inversion of magnetization direc- 
tion (m), which agrees with the expressions proposed 
by Garello et al [see Eqs.(9) and (10) in Ref. [l2|. In 
particular, besides the regular in-plane and perpendicu- 
lar torques captured by Eq. ([1]), additional terms in the 
form oi z X m and m x {z x m) arise. The relative 
magnitude of the different contributions depends on the 
materials parameters. Furthermore, it is interesting to 
notice that such a complex angular dependence of the 
torque is solely determined by the anisotropy in spin- 
relaxation rates (times). By suppressing the anisotropy, 
i.e.; Tz = Txy and 1/t_ = 0, the torque reduces to 
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i+e 



m X y X m 



(7) 



[see Eq.([T])] and the complex angular dependence van- 
ishes. In our model, this anisotropic spin relaxation is 
determined by the D'yakonov-Perel mechanism arising 
from scatterings in the presence of Rashba spin-orbit 
coupling!^ We emphasize here that the above analyti- 
cal results are obtained in a sample of infinite size in the 
strong Rashba coupling regime. 

In the following, we show that the angular dependence 
of the Rashba torque shown in Eq. ([6]) also exists in 
a device of finite size. In addition, we also explain the 
symmetry properties of spin torque at sample edges, as 
shown in Fig. [31 We analyze the angular dependence for 
various m in the xz— plane at three particular locations 
in the device, i.e., in the center (at x,y = 12.5 nm) to 
represent bulk values, and two other locations near the 
edges along y (at a; = 12.5 nm). In what follows, the 
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FIG. 3. (Color online) Angular dependence of spin torque as 
m varies in the xz— plane in a system with L = 25 nm. Panels 
(a) and (c): in-plane torque; panels (b) and (d): out-of-plane 
component. Panels (a) and (b): weak Rashba regime; panels 
(c) and (d): strong Rashba regime. Solid red curves are taken 
at the center of the device. Empty and solid black dots are 
taken at positions near two boundaries along the y direction. 
Insets in panels (a),(b) give a better picture of the solid red 
curves. Solid cyan lines refer to the spin torque for L = cxd. 



description is taken considering Ty and Ti, defined in 
Eq. ©. 

In the weak Rashba regime, when the magnetization 
is along z, the spin density components that contribute 
to the torque show a symmetric profilei^ As m moves 
towards the x axis, the spin density generated perpen- 
dicular to (m, Bfj) points to the z direction. However, 
the faint presence of spin-Hall effect still renders the total 
profile asymmetric, i.e., one edge is more negative than 
the other, in contrast to the case when m = z. Such 
an effect contributes to the angular dependence at the 
edge of the device, as depicted by the open and filled 
black dots in Figl31Ja)(b). In the strong Rashba regime, 
the spin-Hall effect is dominating, producing a more pro- 
nounced angular dependence, as shown by the open and 
filled black dots in FigEl^c) (d) . 

To illustrate the above effects when magnetization is in 
the xz~ plane, wc study the angular dependence in the 
bulk for different device sizes (in Fig. [31 L == 25 nm and 
L ~ CXD. only). In the strong Rashba regime, the spin 
relaxation rate is dominated by D'yakonov-Perel term, 
the angular dependence is pronounced and it does not 
change as the device size increases, which shall eventu- 
ally approach the limit characterized by Eq.((6l). In con- 
trast, in the weak Rashba regime the relaxation rate is 
mostly isotropic, which results in a weak angular depen- 
dence that vanishes as the size increases. These numer- 
ical results are consistent with the argument that the 
angular dependence of spin torque in the bulk is driven 
by the anisotropy in spin relaxation rate. Furthermore, 
in a finite system with isotropic spin relaxation rates, 
oscillations may arise due to edge effects diffusing to- 
wards the center and such a phenomenon is better seen 
in the weak Rashba regime [see inset in Fig. [SJa)]. In 
Fig. [31 the spin torque density in the bulk is plot- 
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FIG. 4. (Color online) Angular dependence as m varies in the 
xz, xy and yz planes in the center of the system for L = 25 nm 
and L — co. Panels (a) and (c): in-plane torque; panels (b) 
and (d): out-of-plane component. Panels (a) and (b) weak 
Rashba regime; panels (c) and (d) strong Rashba regime. Tij 
refers to the spin torque in the ij— plane. The spin torque is 
divided by rUx = cos (j> (ttIz = cos 6) in the xy— {yz—) plane. 



ted when m rotates in the xz—, xy— and yz— planes. 
For the latter two planes, the angular dependence com- 
ing from the cross product in Eq. [T] is removed (divi- 
sion by cos0(cos6')). Therefore, the oscillations of the 
torque magnitude in these three cases can be fitted by 

m 6), which is con- 



Ki + K2 sin^ 0(sin^ 6) + Kg sin"" 



sistcnt with Eq. ([6]). A prominent feature is that the 
oscillations in the a;y-plane vanish as the device size in- 
creases; whereas in the yz-plane, the oscillations persist 
even with an isotropic relaxation rate, which is due to 
F-terni contributing to the angular dependence when m 
is in the j/z-plane, thus enhancing the spin torque ampli- 
tude in the weak Rashba regime [sec Fig. |lUa)(b)]. 

In conclusion, for a finite-size Rashba torque device, we 
have shown that the spin torque is nonvanishing at the 
edges of the sample even when the magnetization and the 
effective Rashba field are parallel, as a result of the com- 
petition between spin-Hall effect and the nonequilibrium 
spin density generated due to anisotropic spin relaxation 
rates. Furthermore, the symmetry and angular depen- 
dence of the spin torque arc in general more complicated 
that the conventional form assumed to date. In our cal- 
culations, the angular dependence is much larger for T\\ 
than for Ti. For a sample of an infinite size, we have ob- 
tained an analytical expression for the spin-orbit torque 
that shows both odd and even components against mag- 
netization inversion and agrees favorably with the expres- 
sion proposed based on experimental results. In a view of 
increasing industrial and academic interests in the field 
of spin-orbit torques, we expect that results presented in 
this Letter shall not only provide a better understanding 
to the key mechanisms behind the experimental observa- 
tions but also shed light on the design of realistic devices. 

We thank I. M. Miron, K. Garello, K. -J. Lee, P. M. 
Haney and M. D. Stiles for valuable discussions. 
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